Coronal Inflows and Giant Polar Plumes 

R. Pinto* ^ R. Grappin* and J. Leorat* 

* Observatoire de Paris, LUTH, CNRS, 92195 Meudon, France 
^Institut d'Astrophysique Spatiale, Universite Paris XI, Or say, France 

Abstract. We present the first results of simulations of giant polar plumes and coronal flows. We use a 2.5D axisymmetric 
MHD numerical model of an isothermal corona and slow solar wind. A plume is generated just above a small magnetic 
bipole embedded in an unipolar flux region which is perturbed by Alfven waves injected from the coronal base. The boundary 
conditions are transparent. The results are compared to those obtained previously with a ID wind model in which plumes are 
generated as a consequence of variations of the heating and flux-tube expansion parameters. 
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INTRODUCTION 

Polar plumes are jet-like overdense features seen over 
the polar regions, mostly during the low solar activ- 
ity phases. These features extend radially, in the direc- 
tion of the magnetic field in the coronal holes. Plumes 
are found to overlie EUV bright points, and these small 
bipoles appear to be undergoing interchange reconnec- 
tion with the unipolar flux concentrations inside coronal 
holes. Outward propagating slow Alfven waves with pe- 
riods between 10 and 15 min and amplitudes of about 
10% — 20% of the plumes' intensity can be observed in 
the first 0.2 Rq d. |[2l[3 | studied the dissipation of 
such waves and its contribution to the coronal heating 
and solar wind acceleration. We focus here in the study 
of the formation and decay of coronal plumes (and asso- 
ciated in/outflows), and report some preliminary results 
of using an MHD numerical model of an isothermal ax- 
isymmetrical solar corona and wind. Alfven waves are 
meant to perturb a small magnetic bipole initially in 
steady-state equilibrium at the bottom of a unipolar coro- 
nal hole. The waves generate unsteady pressure patterns 
which lead to current density accumulation and recon- 
nection. The results are compared to those found by f4l 
for a radially oriented magnetic flux tube which includes 
the dense chromospheric layers and the solar wind. In 
this later case, the wind is perturbed by rapid variations 
of the heating rate in the chromosphere and low corona, 
producing transient in and outflows before the system 
reaches a new steady- state equilibrium. Both cases show 
the formation of an overdense stream which progresses 
outwards in the form of a single wavefront with a phase 
velocity much higher than the bulk stream velocity (of 
the order of the surrounding wind speed). The decay of 
such features is accompanied by transient inflows visible 
mostly in the lower corona. 



Photospheric forcing. Slow photospheric movements 
(null frequency) are expected to shear the coronal mag- 
netic arcades by twisting their footpoints. This "twist- 
ing" may be transmitted upwards and so contribute to 
the magnetic tension/pressure accumulation, for as long 
as the magnetic diffusion remains negligible. Finite fre- 
quency movements, on the other hand, are more prone to 
inject waves {e.g sound and Alfven waves) into the very 
same magnetic structures. The standard numerical ap- 
proach to the modelling of the effects of such movements 
is that of the rigid line-tied forcing of the magnetic loops 
by their footpoints, which are attached to the photosphere 
(both for the null and the finite frequency cases). An im- 
portant consequence of this approximation is that total 
reflection occurs at the foot-points of the loops, which 
then act like resonant cavities; finite-frequency footpoint 
motions induce long-lasting and large amplitude oscilla- 
tions |5| and null frequency motions shear the loop by 
an arbitrarily high amount. The line-tying approximation 
is, nevertheless, one among several possible approxima- 
tions, and an extreme one in some senses. 1 6 1 have shown 
that, at least in the case of a single coronal loop per- 
turbed by "slow" photospheric movements, the line- tying 
approximation may be inappropriate to characterise the 
coronal loop's dynamical response. In what follows, we 
will adopt an alternative approximation which diamet- 
rically opposes the line-tying one: we will consider an 
isothermal MHD corona with a transparent lower bound- 
ary (we'll name it the fully transparent case hereafter, 
for simplicity). We wish to stress at this point that both 
of these are extreme approximations, the physical reality 
laying somewhere in between. The line-tied case is prone 
to severely overestimate the energy input from the pho- 
tosphere into the corona within magnetically connected 
regions, and the fully transparent case, on the other hand, 
underestimates it as the absence of the denser chromo- 
spheric (and photospheric) layers reduces the effective 




FIGURE 1. Initial conditions. Only the northern hemisphere 
is shown, and the radial sub-domain 1 < v/Rq < 3.5. Color- 
scale represents the poloidal Mach number, the dashed contours 
levels \upoi \ /cs = 0.\,\. (with cs ~ 150 km/s) and the white 
lines are magnetic field-lines. The large arrows show where 
Alfven waves will be injected. 



reflection ratios for waves crossing the lower regions. 
12j3\ studied the wave energy leakage through these lay- 
ers which, despite being non-negligible, may be smaller 
than other forms of energy dissipation. The rationale be- 
hind the present model (which uses the fully transparent 
approximation) is to consider a "pessimistic" configura- 
tion in what concerns the destabilisation of coronal struc- 
tures and search for mechanisms which nonetheless are 
able to produce coronal events such as polar plumes, jets 
and inflows. More realistic (less extreme) situations will 
be considered in future work. 

Heating. The ohmic dissipation and consequent heat- 
ing will not be taken into account in our isothermal 
model, but the generation and accumulation of current 
density / can still be quantified. The analysis of the evo- 
lution of / in the domain may give clues to bridge and/or 
compare our results to others from non isothermal mod- 
els. 



GIANT POLAR PLUME 

Numerical model. We use an 2.5 D axisymmetric model 
of the solar corona obeying the following one-fluid 
isothermal MHD equations 
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The magnetic field B decomposes into a time- 
independent external component B^ and an induced 
one b. The solar wind develops into a stable transsonic 
solution in the open field regions. Both the upper and the 
lower boundaries (respectively at r = 14 and 1.01 Rq) 
are transparent. Alfven waves are injected by perturbing 
the corresponding characteristic there. The diffusive 
terms are adapted so that grid scale (A/) fluctuations are 
correctly damped. The kinematic viscosity is defined as 
V = Vo (A//A/o)^, typically with Vq = 2 x 10^"^ cm^ • s~^ 
and 0.01 < (A//A/o)^ < 10. The resistive term is re- 
placed by an implicit filter which dissipates mostly at the 
grid scale and minimises the dissipation of large scales 
fluctuations. Note that actual kinetic dissipation should 
happen at scales much smaller than the grid size, so this 
approximation isn't less realistic than a laplacian term. 
The principals of the numerical model are thoroughly 
discussed in |i9J. The focus here is on circumventing 
the limitations found in O by means of a convenient 
choice of magnetic configuration (by setting B^) and 
Alfven wave injection domain. In particular, we chose 
to perturb regions of the surface not directly connected 
by magnetic field lines, and we do so by starting from 
a magnetic configuration with a null point and separate 
connectivity domains. This way, we can expect wave 
mode coupling to happen due to the highly non-uniform 
phase velocity distributions, the existence of Q = Ca 
interfaces and very large Alfven crossing times near the 
nun point (TOlini. 

Initial conditions. A close-up of the northern polar re- 
gion of the numerical domain is displayed in Fig.[T] This 
particular configuration is comparable to the pseudo- 
streamer structures such as those in 1 12 |. Monochromatic 
Alfven waves will be injected in the zones indicated by 
the large arrows in Fig.[T] The wave period is in the range 
To = 6 — 35 min, the amplitude u^ ^ 49 km/s at the bot- 
tom of the corona and each side will be perturbed with a 
half period phase difference (such that pseudo- streamer 
structure oscillates as a whole). The plasma's j8 is lower 
than 1 all along the lower boundary, and greater than 1 
only over and around the null point. 

Formation of a jet. The Alfvenic wave fronts produce 
a systematic (yet fluctuating) magnetic pressure gradient 
V/?^/2 which push up the denser layers into the j3 > 
1 region, force reconnection and plasma diffusion up 
into the open field region just above the null point. An 
over-dense plasma column forms continuously above the 
null point as the diffused plasma joins the solar wind. 
The jet presents a series of blobs, which correspond 
to density enhancements of amplitude ^ < 0.1 which 
propagate upwards with a phase velocity equal to w// + 
Cs as slow mode wavefronts. (Fig. [2]). The outermost 
edge of the jet progresses also with a velocity equal to 
u// -\-Cs, above the jet's bulk velocity u//. The plasma 





FIGURE 2. Formation of a jet above the bipolar structure in Fig. [T] (turned 30° CCW). The colour-scale represents p corrected 
to compensate for the stratification in the first r ^ 3 Rq of the domain. The time intervals after the beginning of the wave injection 
are, from left to right, 2.5, 4.2, 5.9 and 7.5 h. The amplitude of the injected wave is ^ 49 km/s. The rightmost smaller image 
shows the instant t = 1.5h with tighter colour ranges. 
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FIGURE 3. Profiles of 7^, t/^, p, pUrr^ (top left, top right, 
bottom left and bottom right, respectively) over the white line 
m Fig. [2] Different curves correspond to different instants with 
5t ~ 50 min. 



which does not join the jet falls back along both sides of 
the pseudo- streamer, producing two continuous inflows. 
Fig. |3] shows the evolution of the current density, radial 
velocity, density and mass flux over a cut represented as 
a white line in Fig. [2] A current sheet forms near the null 
point and along the magnetic separatrices as the plasma 
is compressed upwards. Current density peaks around the 
null point. The density contrast between the over-dense 
column and the background grows up to ~ 10 at the base 
of the jet but fades away at higher heights (~ 3 at! Rq 
and ~ 1 at 10 7?0). The plasma's velocity becomes higher 




FIGURE 4. Density n (left) and mass flux nurA (right) as a 
function of time at the bottom of the jet (r = 0.2 ~ 1 .4 R©, the 
same height as for the latitudinal cuts in Fig. [3]). Each curve cor- 
responds to a given frequency coq (correspondent wave periods: 
To = 35, 18, 12, 9, 7, 6 min) for the injected Alfven waves. The 
higher the frequency, the higher the growth rate. The ordinates 
are in log-scale. 



than the surrounding solar wind within the jet, and so 
does the mass flux. 

Growth rate vs. frequency. The Alfven wave driven 
overdense jet may either saturate and reach a quasi- 
steady state or grow continuously for long periods of 
time, depending on the frequency (Oq of the injected 
waves. Fig.|4]shows the growth of the mass flux and den- 
sity within the jet as a function of time for several wave 
frequencies (Oq. There is a smooth transition between a 
stable and an unstable regime at about Tq = ^ = 12 min; 
lower frequencies lead to lower (and eventually null) 
growth rates while higher frequencies lead to ever higher 
growth rates. This behaviour does not seem to break with 
time, at least for as long as the run lasts (Ar ~ 8.5 h). 
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FIGURE 5. Left: peak current density max (7^) /T^. Right: 
JyJ^dV in a volume V which encloses the perturbed region 
(magnetic bipole plus the lower part of the jet). The dotted line 
corresponds to a kinematic viscosity Ivq. 

DISCUSSION AND CONCLUSIONS 

We studied the long-term effects of finite frequency os- 
cillations of the coronal footpoints of a pseudo- streamer 
placed inside a polar coronal hole using an MHD nu- 
merical model of an isothermal corona. We use transpar- 
ent boundaries and show that perturbing such magnetic 
structures at their footpoints with Alfven waves can lead 
to impressive large-scale phenomena. Transverse twist- 
ing movements at the base of the bipolar structure (Fig. 
[T]) lead to a long lasting destabilisation of the system 
and the formation of a dense jet (Fig. [2]) which may 
be interpreted as a polar plume. These movements are 
Alfven waves injected at the base of the bipole, its two 
groups of magnetic arcades being perturbed in phase op- 
position one relative to the other. The effective twisting 
angle of the whole structure remains small (a few de- 
grees) and oscillates with finite frequency, unlike in line- 
tied photospheric shearing models (e.g |13|). A com- 
plex pattern of Alfven wavefronts propagating within 
the lower magnetic arcades non-linearly produces den- 
sity enhancements which are pushed upwards, a part of 
which then diffuses across the magnetic field up into the 
open field region. The diffused plasma then joins the 
overlying polar plume. This diffusion shows a temporal 
modulation with a period equal to that of the injected 
waves (a consequence of exciting both magnetic arcades 
with exactly a half period phase difference; choosing dif- 
ferent phase relations leads to different frequencies). As a 
result, the plume displays a series of blobs which propa- 
gate along its axis as slow mode wavefronts. This result is 
consistent with the slow waves observed by 1 1 1. Current 
accumulates around the null point and along the mag- 
netic separatrices. Fig. [5] highlights the growth of current 
density as a function of time and provides a mesure 
of the production of magnetic small scales. The ohmic 
heating could be defined as Qohm (0 = ^ Sv^^^^^ where 
y encloses the perturbed region 1 14 |. Note that the actual 
values of Qohm depend on the resistive settings and there- 



fore on the local grid-scale, such that Qohm A/~ . Vary- 
ing the kinematic viscosity, on the other hand, as little ef- 
fect on the result. Future work will bring deeper insight 
over the dissipative phenomena. The ohmic heating can- 
not be taken into account here (the model is isothermal), 
but we assume that the simulations of |4|, who use a ID 
non-isothermal model to study the dynamical effects of 
varying the heating rate and flux-tube expansion in the 
lower atmosphere, can represent an open flux tube pass- 
ing close to the magnetic bipole in Fig. [T] The heating 
rate variations at low coronal heights are a proxy to the 
current accumulation and ohmic dissipation. Plume and 
interplume states are also obtained this way, with similar 
predictions for formation and decay time-scales and den- 
sity contrast. On the other hand, the velocity of the wind 
for the plume solution they find is smaller than that for 
the interplume solution. To clear this point, we are work- 
ing on extending the thermal treatment to the 2D case to 
have a unified (and self-consistent) view of these aspects. 
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